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Molecular dynamics simulations are used to study the phase behavior of linear multiblock copolymers with two types of 
monomers, A and B, where the length of the polymer blocks Na and Nb (Na = Nb = N), the number of the blocks ua and 
«b («a = «b = n), and the solvent quality varies. The fraction / of A-type monomers is kept constant and equal to 0.5. Whereas 
at high enough temperatures these macromolecules form coil structures, where each block A or B forms rather individual clusters, 
at low enough temperatures A and B monomers from different blocks can join together forming clusters of A or B monomers. 
The dependence of the formation of these clusters on the varied parameters is discussed in detail, providing a full understanding 
of the phase behavior of linear multiblock copolymers, at least for this symmetrical case. 
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1 Introduction 

The phase separation behavior of block copolymers has been 
the subject of several theoretical and experimental studies.— - — 
In the case of infinitely dilute solutions, it is sufficient to eval- 
uate the behavior of isolated chains. Very interesting is the 
case of multiblock copolymers with blocks composed of either 
A- or B-type monomers for which very successful methods of 
their synthesis exist.- 1,32 Apart from the experimental char- 
acterization, numerical investigations of model chains are the 
most direct approach to understanding the behavior of these 
systems.—"— Interestingly, for the case of multiblock copoly- 
mers there is also close relation to the various toy-models 
(i.e., the HP model— ~^2), which try to mimic the behavior of 
biomacromolecules on the way to understanding complicated 
biological processes, i.e., protein folding, helical structures 43 , 
etc. 

Multiblock copolymers of two chemically different type of 
blocks (A and B) are expected to adopt at high temperatures 
coil structures, where chain conformations are essentially gov- 
erned by repulsive interactions between the different blocks. 
The result is an expansion of the chain dimensions, not only 
with respect to the unperturbed state, but also with respect to 
a homopolymer of the same length under the same thermody- 
namic conditions. 44 Also, one expects that the spherical sym- 
metry of these macromolecules should break, so that the chain 
forms a slightly elongated object, but such phenomena will 
not be discussed in the present manuscript. Here, we focus 
on the microphase separation between the blocks of different 
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type in a multiblock copolymer chain with two types of blocks 
(A and B) and the dependence on its structural parameters in 
terms of their cluster analysis. In this direction, the most inter- 
esting properties of such macromolecules are obtained when 
the solvent quality varies. Under poor solvent conditions (tem- 
peratures lower than the Theta temperature of the chain, which 
in our case will be the same for two different types of blocks 
due to the symmetry of our model) the chain collapses forming 
globular structures, where microphase separation between dif- 
ferent blocks takes place. The different blocks A and B come 
together and form clusters of monomers of type A or B, which 
are microphase separated. In this study, we analyze these phe- 
nomena by studying the formation of these clusters and leave 
the discussion of their size, in conjunction with the size of the 
individual blocks, for a later communication. To this end, we 
have performed large-scale molecular dynamics (MD) simula- 
tions of an off-lattice model of a linear multiblock copolymer 
varying the solvent quality to provide a first approach to the 
understanding of phase separation of such macromolecules. 

The outline of this paper is as follows. Sec. [2] describes 
our model and sketches the analysis needed to characterize 
the size and overall shape (and other properties) of the clus- 
ters considering average properties. Then, Sec. [3]presents our 
numerical results and Sec.|4]summarizes our conclusions. 

2 Model and simulation methods 

In this study we consider the most symmetric case of multi- 
block copolymers, where the fraction of monomers of type 
A is / = 0.5. Then the length of blocks A, Na, and blocks 
B, Nb, is equal as well, and the number of blocks ha with 
monomers of type A and ng with monomers of type B is the 
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same. In this study we denote the total number of blocks as 
n (n = Ha + «b) and the length of the blocks with the capital 
letter N (N = Na = Nb)- Figure Q] shows schematically the 
definition of the above parameters. Then, the total length of 
the multiblock chain is nN. Clearly, this choice facilitates the 
study of this problem, being the most symmetric case. Our 
chains are modelled by the standard bead-spring model,— — 
where all beads interact with a truncated and shifted Lennard- 
Jones (LJ) potential 



U u (r) = 



4eu[(<Wr) 12 




(Pu/r) 



-C ,r<r c 
,r>r c 



(1) 



where r c = 2.5Gzj is the cut-off of the potential, and the 
constant C is defined such that Ujj{? = r c) is continuous at 
this cut-off. Henceforth, units are chosen such that e^j = 1, 
Ou = 1, ks = 1, and m = 1 (mass of the beads) for simplic- 
ity. When we consider two types (A,B) of blocks, we still 
use < = off = off = 1 and efj = eff = 1 but e# = 1/2 
to create an unmixing tendency (Ae = Eab — 1 /2(£aa +£bb))- 
We know that in the case of a binary system with monomers 
at density p = 1 (e.g., a LJ mixture which is a standard sys- 
tem for the study of phase separation), microphase separa- 
tion occurs below a critical temperature T c close to T = 1 .5. 48 
For the multiblock copolymers the average densities are much 
smaller, but since the critical temperature scales proportional 
to the chain length (T c °c %~ r ) t ), we do expect to be able to de- 
tect microphase separation with our model. For this we also 
exploit previous experience with this model where microphase 
separation was studied in the case of bottlebrush polymers 
with two types of side chains. 51 
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Fig. 1 Definition of structural parameters describing our multiblock 
copolymer chains, n is the number of different blocks A and B (in 
this case n = 4) and N is the length of each block. All the blocks, 
irrespective of whether they are of type A or B, have the same length 
N. Then, the total length of the chain is nN. 

The connectivity of the beads along the chain is maintained 
by the finitely extensible nonlinear elastic (FENE) potential 



U F ENE{r) = -\krl\n[\ - (r/r ) 2 ] ,0 < r < r , 



(2) 



where the standard choice of parameters (ro = 1.5 and k = 30) 
was adopted, and UpENE{ r > r o) = °°- 

For the model defined by Eqs. (fl]i and (O the Theta temper- 
ature is known only rather roughly,— namely k 3.0. Being 




Fig. 2 Snapshots of three different multiblock copolymers of the 
same total length nN = 600 at a low temperature T = 1.5 ((a): 
N = 6, (b): N = 15, and (c): N = 60). In case (c) we have the 
formation of two clusters of different blocks which are always phase 
separated. 



interested in T < 0, we have attempted to study the tempera- 
ture range 1.5 < T < 3.0. However, as it is well-known, equi- 
libration of collapsed chains is a rather difficult task and our 
simulation procedure will be discussed below. The tempera- 
ture is controlled by the Langevin thermostat, as done in pre- 
vious work.— ~ 47 i 49 The equation of motion for the coordinates 
{?;(?)} of the beads 



d 2 n 



(3) 



is numerically integrated using the GROMACS packaged In 
Eq. (0 t denotes the time, Ui is the total potential the z'-th bead 
experiences, y is the friction coefficient, and r, (f) the random 
force, y and Ti{t) are related by the fluctuation-dissipation 
relation as follows 



(f i (t)-fj(t')) = 6k B Ty5 ij 8(t-t') 



(4) 



As in previous work,— " 47 ' 49 ' 5 1- — the friction coefficient was 
chosen y = 0.5. For the integration of Eq. Q the leap frog 
algorithm 54 was used with a time step of At = 0.006T, where 
the natural time unit is defined as x = (majj/eu) 1 ^ 2 = E 

As we mentioned earlier, equilibration of collapsed chains 
via MD methods is difficult, and thus we briefly describe here 
our equilibration procedure. First the system was equilibrated 
at T = 3.0 for a time range of 30 x 10 6 x. To gather statistics, 
a sufficient number of statistically independent configurations 
(typically greater than 500) at this temperature was used for 
initial configurations of slow cooling runs, where the temper- 
ature T was lowered in steps of 0.1, and running the system 
at each T for a time which exceeds the relaxation time of our 
chains. The final configuration of each (higher) T was used as 
starting configuration for the next (lower) T. In this way, we 
are able to generate statistically independent configurations of 
clusters. One can see already from snapshot pictures of Fig. [2] 
the different structures that occur at a low temperature T = 1 .5 
for different multiblock chains, where nN is constant, but 
varies. Indeed, rather dense clusters containing several blocks 
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occur, and the conformation depends on the parameters N, n, 
and T (remember that %N « Ae/T). We also note that as the 
length of the chains increases, equilibration of the chains be- 
comes extremely difficult preventing the study of very long 
chains. With our choice of parameters, we expect to observe 
the most interesting effects in the range of simulated chain 
lengths used in this study. In particular, we study in detail the 
chain length nN — 600 and then we discuss the dependence on 
n and N separately, with N of course being the main param- 
eter controlling the incompatibility between A and B, for the 
chosen set of (potential) parameters. 

When the chains collapse at the lower temperatures to form 
microphase separated cluster structures, there is no reason a 
priori for the system to decide which blocks will belong to 
a specific cluster. We actually observe fluctuations where 
a block of type A or B that was part of a cluster of sim- 
ilar monomers escapes from the cluster to become part of 
another neighboring cluster with blocks of the same type of 
monomers. In fact, many such fluctuations would amount to 
influencing the characteristic properties (number N c i of blocks 
a cluster contains, and its linear dimensions) of each cluster. 

For the analysis of the clusters one needs to identify for 
each configuration that is analyzed, which blocks belong to 
which cluster. We have used the standard Stillinger 55 neigh- 
borhood criterion for monomers: if two monomers are less 
than a distance r„ apart, they belong to the same cluster. We 
followed the standard choice r„ = 1.5<7zj and checked that 
qualitatively very similar results were obtained if one chooses 
r„ a bit smaller than this choice (larger values of r n are physi- 
cally hardly significant, since then the particles are too weakly 
bound, due to the rapid fall-off of the LJ potential). The same 
distance r„ is considered also for the analysis of the number 
of contacts. Therefore, a pair of monomers being an abso- 
lute distance less than r„ = 1.5 apart define a "contact". Then, 
the numbers presented in this manuscript denote the average 
number of neighbors per monomer. 

In the case of properties as it is, for example, the number of 
contacts, one just considers the thermal average. However, to 
define properties related to the fluctuating number of clusters 
one should also take into account averages of quantities A that 
depend on these fluctuations, i.e., 

A = ^P(N d )A(N d ). (5) 

As an example, Fig. [3] shows data of the distribution in the 
number of clusters N c i for several choices of n, N. The ex- 
treme choices would be that every block forms a separate clus- 
ter (P(N c f) — 1 for N c i/n = 1, where A and B are always phase 
separated), which is seen for none of the cases presented in 
Fig. [3] (the case n = 10 at temperature T = 3 shows that this 
case is possible with a probability w 0.6), and that all blocks 
of A-type monomers form a cluster with the blocks of B-type 
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Fig. 3 Probability distribution P(N C [) of the number of clusters 
plotted versus N c i /n for combinations of n and N that correspond to 
the snapshots of Fig. [2] Results are shown for two temperatures, 
T = 1.5 and T = 3.0. Lines are guides for the eye. 

forming another cluster (P(N c i ) = 1 for N c i /n — 2 /ri). The lat- 
ter case is taking place for n = 10 at T = 1 .5, where we see a 
single point showing that all blocks of type A are always be- 
longing to the same cluster and all blocks of type B to another 
cluster. The other cases presented in Fig. [3] are intermediate 
ones, where we have a symmetric variation in the number of 
clusters around an average value. 

3 Results and discussion 

The incompatibility between A and B blocks can increase the 
overall chain dimensions of the chains compared to the ho- 
mopolymer chain of the same length under the same thermo- 
dynamic conditions. A and B monomers prefer to be apart, 
increasing in this way the overall size of the macromolecule. 
This can be readily checked if we measure the overall density 
profile of the multiblock copolymers relative to the center of 
mass for various cases of n (n °= 1 /N for nN — const). In Fig.|4] 
we present this dependence on n for nN = 600 = const at a low 
temperature, i.e., T = 1.5. The total chain length (nN — 600) 
is high enough to discuss the different scenarios of microphase 
separation in these macromolecules as the length of the blocks 
varies up to relatively large values of N. Indeed, we can see 
from Fig. [2] that already for N — 60 the formation of single 
clusters takes place. Then we can focus on the variation with 
n and N, with the latter parameter being of course mostly im- 
portant. 

The first quantity we have computed is the overall monomer 
density, mathematically expressed by the following formula 

nN 

P(l?l) = <L S(?- ?c -?,')>, (6) 

;=1 
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Fig. 4 Density profile p(r) plotted versus the radial distance r from 
the center of mass of the chain at temperature T = 1.5. Note the 
pronounced differences occurring for small n, whereas the 
differences for n > 12 are minor. 



where h(x) is the Dirac delta function, r c the position of the 
center of mass of the whole chain, and r, the positions of 
all monomers, irrespective of their type (A or B). The angle 
brackets denote an average over all conformations as usual. 
Figure [4] shows the results for the overall monomer density, 
where n is varied, but nN is kept constant and equal to 600. 
These results refer to a temperature T = 1.5. We can clearly 
see that the homopolymer chain (n = 1) of length (nN — 600) 
obtains always the most compact conformations, if we com- 
pare it with the respective cases of multiblock chains. The case 
of n = 2 exhibits the greatest difference from the homopoly- 
mer case, with the former showing bigger dimensions of the 
chain, due to the formation of a large A-B interface. As n in- 
creases up to n = 12 the multiblock chains gradually exhibit 
more compact structures, but they are always more swollen 
than the homopolymer case. For the range of n between 20 
and 600 shown as an inset in Fig.|4]we see that the differences 
in the density profiles for different multiblock chains are rather 
small. However, for this range of values (« = 20 — 600) we 
observe a monotonic behavior, where the chains obtain less 
compact structures with increasing n. This is clearly due to 
the increase of the unfavorable interactions between A and 
B monomers. Actually, we note here that the increase of n 
should result in breaking of the spherical symmetry, but such 
effects will not be discussed in this manuscript. Here, we 
are only concerned with the phase separation of multiblock 
copolymers under poor solvent conditions. When the density 
profiles of these systems at significantly higher temperatures 
(and closer to the Theta temperature, but still under poor sol- 
vent conditions) are plotted, somewhat smaller differences oc- 
cur and the behavior observed at high N is slightly different. 
These effects will be explained below as we will discuss in 



Fig. 5 Average number of contacts (n con ) per monomer plotted 
versus the number of blocks n for different temperatures 
distinguished with different symbols and line formats. Note the 
logarithmic scale in the x-axis and the different regions that appear 
reflecting the different behavior of multiblock copolymers with 
different n. 



detail the effect of the temperature. 

As n increases, it is natural that the number of contacts be- 
tween A and B monomers increases and then the number of 
contacts A-A and B-B will decrease. Therefore, we are inter- 
ested in the overall number of contacts per monomer irrespec- 
tive of whether a monomer is of type A or B. These results are 
presented in Fig.[5]fbr different temperatures. We observe that 
n con is higher in the case of the homopolymer chain. This con- 
firms our understanding for the homopolymer case, already 
seen in the calculation of the overall profile. The homopoly- 
mer chain is always the most compact globule at this low tem- 
perature T = 1.5 and of course no multiblock chain can obtain 
such a compact structure that the overall contacts are mostly 
favored. Moreover, the results of Fig. [5] show that this is true 
for all temperatures up to T = 3.0. Therefore, the effect of 
incompatibility between A and B is rather weakly present at 
higher temperatures, as expected. These effects would become 
less noticeable as the temperature increases and becomes con- 
siderably higher than the Theta temperature they should be 
unnoticeable. For temperature T = 1.5 and small n the curve 
has a plateau from n = 2 to n = 10. This result is surpris- 
ing, since it shows that the globule at this low temperature has 
bigger dimensions, but the overall number of contacts stays al- 
most constant as the effect of the solvent incompatibility sets 
in strongly. Although there are obvious differences in the den- 
sity profiles of Fig. [4] the multiblock chain for small n rear- 
range in such a way that the number of contacts stays almost 
constant. For values higher than n = 12 the number of con- 
tacts decreases and in the case of very high n(n> 100) we see 
a small increase. In the latter case, we notice that the num- 
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Fig. 6 Average number of clusters N c i divided by the total nunber of 
blocks (n) versus n for different temperatures, as indicated. 



Fig. 7 Number of blocks n\, in each cluster (A or B) divided with 
the number re/2 of A (or B) blocks versus this number (n/2). The 
inset shows the number of blocks in a cluster (A or B) versus n/2. 



ber of monomers per block is small (N = 4- to N = 1). When 
the temperature increases, we see that the plateau for small n 
becomes shorter and for temperatures T > 2.4 has completely 
disappeared. This shows that for T > 2.4 the effect of the tem- 
perature on the behavior of the multiblock chain is minor and 
its behavior is defined by the incompatibility between species 
A and B. We conclude that although the number of contacts 
A-B increases with increasing n, the effect of the temperature 
(when T < 2.4) keeps the number of overall contacts almost 
constant for N > 30. At temperature T > 2.1 the curves show 
a monotonous decrease in the number of contacts and a slight 
increase for high n. At temperature T = 3.0 the variation of n 
in multiblock chains has a very small effect and the difference 
between n — 1 and n = 2 is rather indistinguishable. 

A better understanding of all the above effects can be 
achieved in terms of the analysis of the formed clusters from 
blocks of the same type. In Fig. [6] we have computed the av- 
erage number of clusters with blocks of type A together with 
the different clusters of type B. Of course, when n = 1 we have 
only one cluster with one type of monomers, and for n — 2 we 
have always one cluster with monomers of type A and one 
cluster with monomers of type B and each cluster comprises 
only one block. For small n the ratio N c i jn should be one. 
At temperature T = 1.5, N c i/n stays zero up to n = 12. This 
plateau regime starts at smaller n as the temperature increases, 
in accordance with the results of our previous figures. For 
n > 30 the ratio N c i/n stays almost constant showing that the 
effect of the temperature takes over the effect of the incom- 
patibility between the monomers of different type. At higher 
temperatures, i.e., T = 1.8 and T = 2.1, N c \jn slightly in- 
creases. For even higher temperatures a monotonic behavior 
is observed. The highest variations are observed for high N, 
where the ratio N c i/n decreases rather sharply. 

We can go one step further in order to understand the above 



phenomena. We analyze the number of blocks n/ } of type A (or 
B) in a cluster of type A (or B), respectively. Then, the num- 
ber of blocks A is n/2. Correspondingly the number of blocks 
with monomers of type B is n/2. Due to the symmetry of our 
model we expect the same results for the analysis of clusters 
with blocks of type A and of clusters with blocks of type B. 
Indeed, our results show that this is true confirming the equi- 
libration procedure, which was adopted in this study. In Fig. [7] 
we show such plots where the number of blocks of a cluster 
is plotted versus the total number of A (or B) blocks of the 
multiblock chain. It is now clearly seen the range of N where 
all blocks of the same type belong to the same single cluster. 
This happens when n^/(n/2) = 1. When the temperature is 
rather low (T — 1 .5) this plateau extends up to values n = 20, 
which corresponds to block length ofN = 30. This means that 
the multiblock chain is separated in two clusters, that is one 
of type A and another of type B. This case then corresponds 
to the chain of n = 20, and this explains the results of Figs. [5] 
and [6] When the temperature is lower, the number of blocks 
in the clusters is higher and therefore the curve correspond- 
ing to T = 1.5 shows a smooth gradual decrease. More abrupt 
is the change in higher temperatures where the ratio rij,/ (n/2) 
decreases abruptly. Then the plateau that follows with increas- 
ing n shows the range of n that the enthalpic interactions take 
over. For T > 2.4 the chains exhibit almost the same behav- 
ior with a sharp decrease at small n and a plateau at higher n. 
This shows that at high temperatures and high n hardly ever 
two blocks of the same type come together to form a cluster. 
Then the ratio nt,/(n/2) w 1/ {n/2). However, at temperature 
T = 1.5 this ratio has the value n\j/(n/2) w 0.2 showing that 
the formation of clusters with more than two blocks is possi- 
ble. The inset shows the above results in a different represen- 
tation. Then the linear curve (y = ax) denotes the regime that 
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Fig. 8 Standard deviation in the number of blocks in each cluster (A 
or B) divided with the number n/2 of A (or B) blocks that the 
multiblock chain contains versus this number (n/2). The inset 
shows the standard deviation in the number of blocks in a cluster (A 
or B) versus n/2. 



the ratio rib/ (n/2) = 1. It is also seen that this regime becomes 
smaller as the temperature is increased. Then, at low temper- 
atures, the incompatibility between the monomers of different 
type becomes significant. This is not seen at high temperatures 
close to Theta, where the effect of the solvent quality and the 
enthalpic interactions rather cancel each other. 

The number of blocks contained in a cluster is not the same, 
but it is a quantity that fluctuates, as we have already seen from 
Fig. [3] Then Fig. [8] shows these fluctuations of the number of 
blocks in a cluster. For small values of n (or high N) such 
fluctuations should not exist, since there are only two clusters, 
each one with blocks of monomers of different type. For in- 
termediate values of n the fluctuations in the number of blocks 
in the clusters become rather high. The blocks of type A try 
to rearrange in such way that they form a single cluster, but 
the number n is high enough to hinder such effect. For n in 
the regime of strong fluctuations the formation of single clus- 
ters containing all blocks of the same type is possible. For 
higher n the curves show a monotonic smooth behavior and 
the enthalpic interactions take over. It is also noticeable that 
the most pronounced fluctuations are observed in lower tem- 
peratures, where the effect of the incompatibility between dif- 
ferent monomer types and the effect of the solvent incompat- 
ibility compete. At T = 1.5 the fluctuations persist over high 
values of n. This is clearly seen from the inset of Fig. [8] Also, 
as N becomes smaller, the fluctuation in the number of blocks 
in a cluster is higher. 

Interesting is the dependence on the temperature for this 
particular case. From Fig. [9] we can clearly distinguish three 
different cases. For small values of n (2 < (n/2) < 6) the 
formation of a single cluster is possible up to temperatures 
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Fig. 9 Number of blocks ny in each cluster (A or B) divided with 
the number n/2 of A (or B) blocks that the multiblock chain 
contains versus temperature. 



T = 1.8. In this case the blocks of the same type belong al- 
ways to the same cluster. This behavior is weakened as n in- 
creases, or correspondingly N decreases. The second case is 
that for 10 < (n/2) < 30, where the formation of clusters con- 
taining all blocks of the same type is possible. Of course, this 
is more probable when N is higher and the temperature is low. 
As the temperature is increased, we obtain a strong decrease 
in the number of blocks per cluster and at high temperatures 
T > 2.3 we obtain a plateau where n^/(n/2) ss 2/n. When n 
is high the dependence on the temperature is rather small and 
the behavior of the chain is defined by the incompatibility be- 
tween A and B monomers. Of course for the latter case the 
formation of single clusters with all blocks of the same type is 
impossible. 

Figure [10] shows the dependence of the fluctuation in the 
number of blocks per cluster. The stronger fluctuations are 
observed at lower temperatures (depending on n) in the range 
between T = 1.7 and T = 2.4. As n increases this peak of the 
higher fluctuations is shifted to the left at lower temperatures. 
The results of Figs.|9]and[T0]show that at low temperatures the 
formation of single clusters is possible, but the same time the 
fluctuations in the number of blocks in the cluster are rather 
high. 

For block copolymers we know that the incompatibility be- 
tween different monomers is expressed by the term %N. Here, 
we discuss the effect of the block length N. We choose three 
different chain lengths, which correspond to low N, interme- 
diate, and high enough, where always the blocks of the same 
type form a single cluster. Then we look for the dependence 
on n. The results of Fig. QT| confirm the above argument that 
the phase separation depends strongly on N, as expected, and 
rather weakly on n for a collapsed globule. Of course the in- 
fluence of the surface energy (°= R 2 , with R being the radius 
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Fig. 10 Standard deviation in the number of blocks n/, in each 
cluster (A or B) divided with the number n/2 of A (or B) blocks that 
the multiblock chain contains versus temperature. 
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Fig. 11 Number of blocks n/, per cluster divided with the number of 
blocks of type A (or B) versus this number of blocks n/2. Results 
concern different block lengths N and temperatures T. 



of the almost spherical collapsed object) is a combination of 
the above two parameters, n has such effect with N having the 
additional effect of controlling the incompatibility between A 
and B monomers. When N is high (N = 80 or 160), the blocks 
are joined together into a single cluster with monomers of the 
same type at temperature T = 1.5. Therefore, for this case 
no dependence on n is observed. However, when n becomes 
larger than some value, the formation of more than one cluster 
takes place. Then the overall length of the multiblock chain 
is very long and difficult to study with simulation techniques. 
Thus, such effect is not discussed in this manuscript and is 
not systematically studied. For intermediate N (N = 20) and 
small n the formation of a single cluster is still possible, but 
this length is not high enough to ensure strong incompatibility 
between the blocks with monomers of different type. Some- 
what richer behavior is expected as it has already been dis- 
cussed, with strong fluctuations in the number of blocks per 
cluster. For N = 5 the formation of a single cluster is not 
seen and the number of blocks per cluster and per block de- 
creases reaching a plateau value at high n. For T = 3.0 we see 
that that the different block length does not play any role. The 
three different block lengths have collapsed into a single curve 
(%N <x Ae/T), This is the indication that the behavior of our 
multiblock chain depends strongly on the effect of the temper- 
ature. At high temperature also the curve reaches a plateau, 
that is for (n/2) > 25. 

Figure [12] shows the dependence of the number of blocks 
per cluster on N. For (n/2) = 4 and at low temperature we 
have the formation of a single cluster for N > 20. As the tem- 
perature increases this happens for higher values of N. At 
T = 2.2 the curve takes gradually a constant value. For tem- 
perature T = 3.0, which corresponds roughly to the Theta tem- 
perature, there is no dependence on N, as it is expected. This 



is also observed for n = 20 and n = 100. Comparing the re- 
sults for different n at the low-temperature regime, we see that 
the phase separation depends rather strongly on N and rather 
weakly on n. We should note here that all phase changes oc- 
cur gradually due to the presence of thermal fluctuations in a 
collapsed object with a high ratio of surface to bulk energy. 
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Fig. 12 Same as Fig.QTJ but the dependence on the block length N 
is shown. The curves correspond to different numbers of blocks n 
and temperatures T . 



4 Conclusions 

In this manuscript we have presented a detailed study of the 
phase behavior of multiblock copolymers varying the number 
of blocks n, the length of the blocks N, and the solvent qual- 
ity by variation of the temperature. We have described the 
different types of phase separation appearing in linear multi- 
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block copolymers and we have discussed in detail the effect of 
these parameters. Although the temperature (solvent incom- 
patibility) and the block length N, as in the bulk, control the 
phase separation between blocks of different type (for given 
Ae), the effect of the number of blocks n is also important 
due to the increase of the surface free energy. In the range of 
values that have been discussed here, very reach behavior is 
observed for intermediate values of «, N and T, where the fluc- 
tuations in the number of blocks per cluster are pronounced. 
In this case the interplay between enthalpic and entropic ef- 
fects compete resulting in different types of microphase sep- 
aration in the collapsed states. The behavior of multiblock 
copolymers can be parallelized with that of various biological 
macromolecules which are formed by periodically repeated 
structural units ("monomers") along their chain. Our study 
provides a first step towards the understanding of the phase 
separation of such complex systems. 
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